function e = strakos(n,lambda1, kappa, rho)
e=zeros(n,1);
e(1)=lambda1;
for i=2:n,
    e(i)=lambda1+(i-1)/(n-1)*(kappa*lambda1-lambda1)*rho^(n-i);
end